Journal  of  Power  Sources  185  (2008)  241-247 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


Fractal  model  for  prediction  of  effective  thermal  conductivity  of  gas 
diffusion  layer  in  proton  exchange  membrane  fuel  cell 

Ying  Shia’b,  Jinsheng  Xiaobc’*,  Shuhai  Quana,  Mu  Panb,  Runzhang  Yuan'3’* 1 * 

a  School  of  Automation,  Wuhan  University  of  Technology,  Hubei  430070,  China 
b  State  Key  Laboratory  of  Advanced  Technology  for  Materials  Synthesis  and  Processing, 

Wuhan  University  of  Technology,  Hubei  430070,  China 

c  School  of  Automotive  Engineering,  Wuhan  University  of  Technology,  Hubei  430070,  China 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  2  July  2008 
Accepted  9  July  2008 
Available  online  17  July  2008 


Keywords: 

Proton  exchange  membrane  fuel  cell 

Gas  diffusion  layer 

Effective  thermal  conductivity 

Prediction 


The  process  of  heat  transfer  within  porous  media  is  usually  considered  as  a  transport  through  large 
numbers  of  straight  channels  with  uniform  pore  sizes.  For  the  prediction  of  effective  thermal  conductivity 
of  gas  diffusion  layer  (GDL),  morphological  properties  such  as  the  tortuosity  of  channels  and  pore-size 
distribution  of  this  porous  layer  should  be  considered.  Thus  in  this  article,  novel  parallel  and  series-parallel 
prediction  models  of  effective  thermal  conductivity  for  the  GDL  in  proton  exchange  membrane  fuel  cell 
(PEMFC)  have  been  derived  by  fractal  theoretical  characterization  of  the  real  microstructure  of  GDL.  The 
prediction  of  fractal  parallel  model  for  carbon  paper,  a  basal  material  of  the  GDL,  is  in  good  agreement 
with  the  reference  value  supplied  by  Toray  Inc.  The  prediction  results  from  the  proposed  models  are  also 
reasonable  because  they  are  distributed  between  the  upper  and  lower  bounds.  Parametric  effect  has  been 
investigated  by  using  the  presented  models  in  dimensionless  formalism.  It  can  be  concluded  that  dimen¬ 
sionless  effective  thermal  conductivity  (k(,ff)  has  a  positive  correlation  with  effective  porosity  (e)  or  the 
pore-area  fractal  dimension  (Dp)  when  fcs/kg  <  1 ;  whereas  it  has  a  negative  correlation  with  e  or  Dp  when 
ks/kg>l  and  with  tortuous  fractal  dimension  (Dt)  whether  ks/kg<  1  or  not.  Furthermore,  these  fractal 
models  have  been  modified  by  considering  the  effect  of  polytetrafluoroethylene  (PTFE)  incorporated  into 
the  pore  spaces  of  carbon  paper,  and  the  corresponding  model  prediction  shows  that  there  is  an  increase 
in  the  effective  thermal  conductivity  due  to  the  filling  of  PTFE  that  has  high  thermal  conductivity. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Thermal  conductivity  is  generally  considered  as  a  significant 
material  parameter  in  the  heat  transfer  within  a  gas  diffusion 
layer  (GDL),  which  is  modeled  as  a  uniform  porous  medium, 
and  is  inevitably  used  in  the  thermal  simulation  of  a  proton 
exchange  membrane  fuel  cell  (PEMFC).  In  this  article,  the  pro¬ 
cess  of  heat  transfer  within  porous  GDL  has  been  investigated  and 
suitable  prediction  models  of  effective  thermal  conductivity  have 
been  suggested  by  fractal  theoretical  characterization  of  the  real 
microstructure  of  GDL. 

The  complexity  of  heat  transfer  in  inhomogeneous  porous  media 
is  due  to  both  solid  and  fluid  phases  and  random  pore  morphology. 
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The  process  of  heat  transfer  follows  different  patterns  for  varied 
pore  structures.  In  general,  heat  transfer  is  controlled  by  the  mech¬ 
anisms  of  convection,  radiation,  and  conduction. 

Convective  heat  transfer  occurs  when  fluid  flows  within  pores. 
Generally,  the  effect  of  convection  is  more  apparent  in  the  case  of 
large  pore  sizes,  whereas  it  can  be  neglected  for  small  pore  sizes 
(<100  p.m)  at  lower  temperatures  (<373 1<)  [1]  because  of  lack  of 
intensive  fluid-circulation  within  the  pores. 

Radiative  heat  transfer  occurs  through  the  adsorptive  or  radia¬ 
tive  heat  emission  of  pore  walls.  From  a  study  on  coal  chars,  it  has 
been  found  that  the  mechanism  of  radiation  has  a  significant  effect 
on  heat  transfer  through  large  pore  sizes  (>10  pan)  at  high  temper¬ 
atures  (>1000  K)  [2],  whereas  for  most  carbonaceous  materials,  the 
effect  of  radiative  heat  transfer  can  be  neglected  for  temperatures 
below  1000  K. 

For  the  GDL  of  PEMFCs,  because  the  working  temperature  is 
lower  than  373 1<  and  because  the  pore  sizes  of  GDL  are  smaller 
than  100  pan  [3],  neither  convection  nor  radiation  is  significantly 
involved.  Thus,  heat  transfer  within  GDL  is  controlled  only  by  the 
mechanism  of  conduction.  Although  the  process  comprises  only 
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heat  conduction,  it  is  still  complex  because  it  involves  the  following  on  porous  media  with  idealized  structures,  there  are  two  types  of 
three  processes:  theoretical  models,  effective  medium  model  and  layer  model. 


(1)  Heat  transport  within  carbon  fiber; 

(2)  Heat  transport  within  fluid  in  pores; 

(3)  Heat  transport  between  connected  carbon  fibers. 

The  dominant  mode  of  transfer  is  determined  by  the  contribu¬ 
tion  of  the  solid  and  fluid  phases  to  thermal  conductivity.  If  the 
thermal  conductivity  of  the  fluid  is  greater  than  that  of  the  carbon 
fibers,  the  heat  transfer  through  fluid  becomes  dominant.  The  com¬ 
plexity  of  the  heat-transfer  process  is  attributed  to  various  factors, 
such  as  material  properties  of  the  solid  phase;  shape,  size,  and  size- 
distribution  of  pores;  and  type  of  component,  state,  and  properties 
of  fluids.  In  addition,  pressure  and  temperature  also  have  signifi¬ 
cant  effects  on  heat  transfer.  Therefore,  the  exact  determination  of 
the  effective  thermal  conductivity  is  highly  difficult. 

2.  Prediction  methods  for  effective  thermal  conductivity  of 
porous  media 

In  general,  thermal  conductivity  can  be  predicted  by  empiri¬ 
cal  formulas,  numerical  simulations,  or  theoretical  models.  Table  1 
shows  some  specific  values  of  thermal  conductivity,  which  are 
derived  from  experimental  measurements  [4], 

2.1.  Empirical  formulas 

Most  of  the  pure  empirical  or  semiempirical  formulas  with  sev¬ 
eral  empirical  constants  are  derived  by  polynomial  fitting  methods 
[2,5,6].  Generally,  these  empirical  constants  do  not  indicate  any 
specific  physical  meanings,  and  their  values  derived  from  different 
researches  lack  uniformity. 

2.2.  Numerical  simulations 

Numerical  simulations  can  be  used  for  the  research  of  het¬ 
erogeneous  media  with  discrete  hierarchies.  For  the  study  of  a 
porous  medium,  it  is  necessary  to  build  a  geometrical  model  of 
the  porous  medium  and  specify  one  of  the  numerical  methods, 
such  as  finite  element  method  [7],  finite  difference  method  [8],  or 
boundary  element  method  [9].  Subsequently,  through  simulation, 
thermal  conductivity  can  be  predicted  and  parametric  effect  can 
be  studied.  However,  because  the  geometrical  model  varies  with 
the  microstructure  of  porous  media,  it  always  needs  rebuilding  for 
different  research  objects;  moreover,  numerical  simulation  is  time- 
consuming,  and  this  prediction  method  is  unsuitable  for  fast  and 
convenient  prediction  of  thermal  conductivity. 

2.3.  Theoretical  models 

Theoretical  model  prediction  of  thermal  conductivity  is  carried 
out  by  using  theoretical  expressions.  For  the  theoretical  researches 


Table  1 

Various  specific  values  of  thermal  conductivity 


Material  Thermal  conductivity  Material  Thermal  conductivity 

(Wm-'K-1)  (Wm-'K-1) 


Silver  4.186  x  102 

Aluminum  2.093  x  102 

Quartz  8.392 

Sandstone  3.767 

Clay  0.837-1.256 

Water  0.461 


Hydrogen 

Oxygen 

Air 

Benzene 

Petroleum 


0.167 

0.025 

0.026 

0.159 

0.147 

0.502-l.C 


2.3.1.  Effective  medium  model 

Typical  effective  medium  models  are  sphere,  ellipsoid,  cube,  and 
tube  models. 

Hashin  and  Shtrikman  [10]  have  applied  the  sphere  model,  along 
with  the  variation  principle,  to  find  the  upper  and  lower  bounds 
of  the  effective  magnetic  conductivity  for  a  multiphase  mixture 
on  the  basis  of  uniform  and  isotropic  assumptions.  These  bounds 
are  also  applicable  to  the  thermal  conductivity.  Furthermore,  the 
lower  bound  /qow  and  upper  bound  kUpper  can  be  determined  by 
the  following  equations: 

kiow  =  ki  +  1/(j£2_kij  +  ei/3ki  (1) 

kupper  =  k2+1/(ki_kJ;  +  £2/3k2  (2) 

where  k2>kf,  they  are  the  widely  accepted  bounds  for  uniform 
and  isotropic  materials  and  are  dependent  only  on  the  volumetric 
fraction. 

Grant  and  West  [  11  ]  have  adopted  a  three-dimensional  tube  unit 
to  emulate  a  two-phase  system  and  have  obtained  the  effective 
thermal  conductivity  ke{{  as 

keff=lxmkm  +  (l-Xm)ks  (3) 

where  km  and  Xm  are  the  thermal  conductivity  and  the  volumetric 
fraction  of  tubular  inclusions,  respectively,  and  ks  is  the  thermal 
conductivity  of  the  second  phase. 

Hsu  etal.  [12]  have  used  a  three-dimensional  lumped-parameter 
model  for  the  derivation  of  the  effective  thermal  conductivity  of 
spatially  periodic  porous  media  and  have  obtained  the  following 
result: 


keff 

kf 


1  -  7a  -  2 VaYc  +  2 yIyc  + 


vl  -  yIyI 

1  -  Ka  +  Ya^ 


,  2(yayc  -  YaYc) 

1  —  YaYc  +  YaYc^- 


(4) 


where  kf  is  thermal  conductivity  of  the  fluid  phase,  ya  =  a//,  yc  =  c/a, 
/  and  a  are  the  geometrical  lengths,  c  is  the  contact  length,  A.  repre¬ 
sents  fluid/solid  thermal  conductivity  ratio. 

Every  parameter  in  these  effective  medium  models  has  a  clear 
physical  meaning,  but  it  is  difficult  to  obtain  some  of  the  model 
parameters;  moreover,  the  above  models,  not  the  general  ones,  can 
be  only  applied  to  specific  porous  media. 


2.3.2.  Layer  model 

The  parallel  or  series  layer  models  are  widely  used  for  prediction 
of  the  effective  thermal  conductivity  of  a  two-phase  system.  The 
effective  thermal  conductivity  kp  can  be  calculated  by  the  parallel 
model  using  the  following  equation: 


kp  =  kfEf  +k2E2 


(5) 


where  s\  and  s2  are  the  volumetric  fractions,  respectively;  and  kf 
and  k2  are  the  effective  thermal  conductivities  of  the  two  phases, 
respectively. 

Using  the  series  model,  effective  thermal  conductivity  ks  can  be 
obtained  as 


1 


~  =  ir  +  ir 


(6) 


Both  parallel  and  series  models  are  usually  used  for  the  deter¬ 
mination  of  the  upper  and  lower  bounds  of  the  effective  thermal 
conductivity. 
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For  porous  media,  numerous  capillary  channels  both  parallel 
and  perpendicular  to  the  flow  direction  (Fig.  1 )  exist;  therefore,  the 
effective  thermal  conductivity  can  be  determined  by  the  series- 
parallel  model: 


keff  ((i-5)/kp)+tt/fcj  (7) 

where  f  is  the  ratio  of  the  number  of  perpendicular  channels  to  the 
total  number  of  channels. 

Usually,  ideal  geometrical  models  are  used  for  inferential  rea¬ 
soning  formulas.  Among  these  models,  the  straight  capillary  model 
with  uniform  pore  sizes  is  widely  used.  Because  these  geometrical 
models  are  regularly  arranged,  they  are  inapplicable  to  real  porous 
media  with  disordered  microstructural  characteristic. 

Thompson  et  al.  [13]  have  found  that  real  porous  media  have 
a  fractal  characteristic.  Several  studies  have  focused  on  the  appli¬ 
cation  of  the  fractal  method  to  the  prediction  of  effective  thermal 
conductivity  [14-19].  It  can  be  accomplished  by  two  methods.  One 
method  is  the  application  of  fractal  geometrical  models  such  as 
random  Sierpinski  carpets  in  the  simulation  of  thermal  equations 
[14,15],  In  contrast  with  the  normal  volume-average  models,  the 
fractal  models  are  more  realistic  because  they  consider  the  disor¬ 
dered  characteristic,  but  the  numerical  simulation  needed  in  this 
method  is  time-consuming.  The  second  method  is  the  application 
of  fractal  theoretical  models,  along  with  the  effective  medium  and 
layer-analog  methods  [15-19].  This  method  is  suitable  for  both  ana¬ 
lytical  solution  and  fast  prediction.  Therefore,  this  method  has  been 
adopted  here. 

A  previous  study  by  the  authors  has  shown  that  microstructures 
and  pore-size  distribution  of  the  GDL  in  PEMFC  have  a  fractal  char¬ 
acteristic  [3],  In  the  following  section,  the  analytical  solution  model 
of  the  effective  thermal  conductivity  is  derived  on  the  basis  of  the 
fractal  description  of  GDL.  Subsequently,  the  effects  of  microstruc¬ 
tural  parameters  on  effective  thermal  conductivity  are  investigated. 


3.  Fractal  model  of  effective  thermal  conductivity  of  GDL  in 
PEMFC 


3.1.  Fractal  parallel  model 


Experimental  investigations  on  permeation  in  porous  media 
have  shown  that  the  channels  through  which  liquids  permeate  have 
fractal  characteristics.  Similarly,  the  heat-transfer  routes  within 
porous  media  may  also  have  fractal  characteristics  and  can  be  rep¬ 
resented  as  random  Koch  curves.  Porous  carbon  paper,  made  from 
graphite  fibers,  is  used  for  the  preparation  of  the  GDL  in  PEMFC.  The 
heat  transfer  within  GDL  is  assumed  to  be  similar  to  that  within 
the  tortuous  fractal  parallel  channels  with  different  pore  sizes  and 
is  depicted  in  Fig.  2. 

The  length  L(  k )  of  the  capillary  pathway  is  related  to  the  capillary 
size  k  (i.e.  the  pore  diameter)  by  the  following  fractal  relationship 
[20]: 

where  L0  is  the  representative  or  linear  length  of  these  capillary 
pathways  towards  the  flowing  direction,  and  Dt  is  the  tortuous 
fractal  dimension. 

The  tortuosity  of  the  capillary  pathway  r  can  be  obtained  by 


(9) 


Another  characteristic  of  porous  media  is  that  the  cumulative 
pore  population  N  in  a  unit  cross-section  may  be  mathematically 
expressed  as  follows  [20]: 

N(L>k)=  (5||)DP  (10) 


where  Dp  is  the  pore-area  fractal  dimension,  and  k  and  A,max  are  the 
pore  size  and  the  maximum  pore  size  of  porous  media,  respectively. 
The  first  derivative  of  Eq.  (10)  with  respect  to  k  is 

-dN  =  DpAS,axA-(DP+1)  dk  (11) 


The  total  area  of  the  unit  cell  A  can  be  described  as  follows  [21  ] : 


7rD,,/,(nilx  r 
=  48(2  —  Dp)  [ 


(12) 


On  the  basis  of  the  above  fractal  descriptions  of  the  GDL,  the  heat 
transfer  through  parallel  channels  is  studied,  and  than  the  proposed 
model  is  modified  by  the  series-parallel  method.  A  schematic  rep¬ 
resentation  of  a  porous  medium  with  parallel  channels  is  shown 
in  Fig.  3,  where  s  is  the  porosity  and  A  is  the  cross-sectional  area 
perpendicular  to  the  direction  of  heat  flow;  then,  the  effective  area 
of  heat  transfer  for  the  solid  phase  is  sA.  The  heat  transfer  within  a 
porous  GDL  can  be  simplified  as  a  transport  through  parallel  chan¬ 
nels  of  the  solid  and  gas  phases  as  shown  in  Fig.  3,  where  the  parallel 
gas  channels  have  fractal  structure,  and  Rg  and  Rs  are  the  ther¬ 
mal  resistances  of  the  gas  and  solid  phase,  respectively,  k%  and  ks 
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the  GDL. 


are  their  respective  thermal  conductivities,  L0  indicates  the  linear 
lengths  of  the  channels  for  the  transfer  of  heat  flux  Q,  and  L(X) 
denotes  the  actual  length  of  the  tortuous  channel. 

First,  a  fractal  parallel  model  of  effective  thermal  conductivity 
for  the  GDL  is  deduced  as  follows. 

The  thermal  resistance  of  a  single  channel,  R(X),  can  be 
expressed  as 


I(A.)  L(X)  4L(X) 

1  J  A(X)kg  (nX2/4)kg  nX2kg 


(13) 


Furthermore,  the  number  of  channels  with  pore  size  X  is  -d N;  their 
parallel  thermal  resistance  RN(X )  can  thus  be  expressed  as 


Rn(X ) = 


RW 

-dN 


AL(X) 

7zX2kgDpXty3xX-(DP+'1')  dX 


(14) 


keff,p  — 


JtkgDfX*£[  1  -  (A.min/^max)Dt'Dp+1  ] 
4Lgt~1(Dt  —  Dp  +  1) 


4(2  -  Dp)s 

x  7rDpX'iax{i  -  (xmm/xmaxy-un 

kg(2  -  Dp)sX^[\  -  (Amin/Amax) 
~  L^-\Dt-Dp  +  \)[\-(Xmin/XmB 


+  (1  -  £)ks 

Dt— Dp+1 1 

1X^2-DP|  +  (1  -  £)ks 


(20) 


This  is  the  fractal  model  of  effective  thermal  conductivity  for  the 
GDL  based  on  parallel  theory. 

When  the  channels  are  straight,  i.e.  Dt  =  1,  substituting  this  value 
into  Eq.  (20),  the  following  relation  is  obtained: 


keff,P  =  kg£  +  (l  -e)fcs  (21) 

It  is  the  general  formula  for  the  calculation  of  effective  thermal 
conductivity  for  a  two-phase  parallel  system. 

The  comparison  of  Eq.  (20)  with  Eq.  (21)  indicates  that  their 
second  terms  are  the  same,  whereas  the  first  ones  are  not.  Eq.  (20) 
shows  that  the  microstructural  parameters  affect  only  the  thermal 
conduction  of  the  gas  phase  because  the  fractal  characteristics  of 
the  pores  that  have  an  effect  only  on  gas  phase  are  considered. 

From  Eq.  (20),  it  also  can  be  observed  that  the  parameters  have 
clear  physical  meanings  and  there  are  no  empirical  constants.  This 
general  theoretical  model  can  be  used  both  for  prediction  of  effec¬ 
tive  thermal  conductivity  of  two-phase  porous  systems  and  for  the 
study  of  the  effect  of  porosity  and  other  microstructural  parameters 
(such  as  tortuous  fractal  dimension  Dt  and  pore-area  fractal  dimen¬ 
sion  Dp).  From  the  viewpoint  of  characterization  of  microstructure, 
this  fractal  model  is  closer  to  reality  than  the  common  models  based 
on  volume  average. 


3.2.  Modification  of  the  fractal  parallel  model 


Substituting  the  definition  formula  of  L(X)  into  Eq.  (14),  the  follow¬ 
ing  relation  is  obtained: 

4L^tX1~°t  4Lqc 

=  7TkgDpA.KfaxA.i-Dp  dX  =  JtkgDpXfyxXDt-Dp  dX  (15) 

According  to  the  parallel  principle  of  electric  resistance,  the  thermal 
resistance  of  the  gas  phase,  Rg,  can  be  calculated  by  the  expression 


Rg  1  /RN(X)  J'2X(xkgDpX±xX°t-»r/4L»‘)  dX 

4Lpt(Dt  —  Dp  +  1) 

“  7tkeDpX%£  [1  -  (A.min/A,niax)Dt_Dp+  ‘  ] 


(16) 


The  thermal  resistance  of  the  solid  phase,  Rs,  is  as  follows: 


Rs 


Lo 

(1  -  e)Aks 


(17) 


Thus,  the  total  resistance  of  the  porous  medium,  R,  can  be 
obtained  as  follows: 


1  =  L0 

(1/Rg)  +  (1/J?s)  Ake  ffiP 


(18) 


where  keff  p  is  the  parallel  effective  thermal  conductivity  of  the 
porous  GDL,  and  from  Eq.  (18),  the  following  expression  is  derived: 


nkgDpX%+'  [1  -  (A.minAmax)Dt~Dp+1  ] 

=  4ALq'-  '  (Dt  -  Dp  +  1 )  H 


Substituting  Eq.  (12)  into  Eq.  (19),  ke{{  p  is  obtained: 


Because  there  are  channels  both  parallel  and  perpendicular  to 
the  direction  of  heat  flow  within  GDL  in  PEMFCs,  both  their  effects 
should  be  considered.  For  the  ultrathin  GDL  with  a  thickness  of 
1.9  x  10-4  m,  there  may  be  a  few  channels  perpendicular  to  the 
direction  of  heat  transfer.  The  series  theory  can  be  applied  to  modify 
the  fractal  parallel  model  proposed  above  and  to  obtain  the  effective 
thermal  conductivity  of  the  gas  and  solid  phases  in  series: 


keff,s  = 


1 

£/kg  +  (l  -  s)/ks 


(22) 


Using  the  series-parallel  model,  the  effective  thermal  conductivity 
of  the  GDL  can  be  expressed  as 


keff  = 


1 

((l-?)/keff,p)  +  (?/keff,s) 


(23) 


where  £  is  the  ratio  of  the  number  of  perpendicular  channels  to  the 
total  number  of  channels,  with  values  ranging  from  0  to  1. 

Eqs.  (20),  (22),  and  (23)  can  be  rewritten  in  terms  of  dimension¬ 
less  variables,  by  dividing  both  sides  of  these  equations  by  kg: 

keff.p  (2  - Dp)eX%-J [1  -  (AlinAmaX)Dt~Dp+1  ] 
eff’p  kg  Lgt-^Dt  -  Dp  +  1)[1  -  (A.minAmax)2_Dp] 

+  (l-£)^  (24) 

. ,  _  keff,s  _ _ 1 _ 

eff’s  kg  £  +  (l  -e)/(ks/kg) 

j./  _  keff  _ _ 1 _ 

eff“  kg  “((H)^)^ 


(19) 


(26) 
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In  the  above  derivations,  the  GDL  is  considered  as  a  two-phase 
system  and  a  general  fractal  model  has  been  deduced.  Because 
carbon  paper,  a  basal  material  of  GDL,  is  always  treated  with  polyte- 
trafluoroethylene  (PTFE)  before  being  applied  to  PEMFC,  the  effect 
of  the  PTFE  incorporated  into  the  pore  spaces  of  carbon  paper 
should  be  considered  as  follows: 

,  ks(2  -  Dp )£^t-i  [ !  _  ( W  W)D‘-De+1  ] 

eff’P  Lq'-1  (Dt  -  Dp  +  1 )[  1  -  (Amin  Amax)2~Dp  ] 

+  0.22fcsl  +  (0.78  -  s)kS2  (27) 


keff’s  s/ks  +  0.22 /ksl  +  (0.78  -  s)/ks2  (28) 

keff=((l-|)/feeff,p)  +  (?/fceff,s)  (29) 

where  k5]  and  ks2  are  the  thermal  conductivities  of  the  carbon  fiber 
and  PTFE,  respectively.  The  modified  equations  are  suitable  for  the 
prediction  of  the  effective  thermal  conductivity  of  a  GDL. 

Eqs.  (23 ),  (26),  and  (29)  are  all  functions  of  the  fractal  dimensions 
Dp,  Dt,  and  other  microstructural  parameters,  among  which  there 
is  only  one  empirical  constant  with  a  clear  physical  meaning. 

4.  Results  and  discussion 


4.1.  Fractal  prediction  of  effective  thermal  conductivities 

The  effective  thermal  conductivities  ke{{  can  be  predicted  by  Eq. 
(23)  or  (29)  when  air  is  passed  through  the  samples  (two  sam¬ 
ples  are  used:  sample  ‘a’  is  the  TGP-H-060  carbon  paper;  sample  ‘b’ 
is  the  TGP-H-060  carbon  paper  treated  with  PTFE).  Subsequently, 
a  comparison  is  made  between  different  theoretical  models.  The 
parameters  of  GDL  used  in  the  study  are  listed  in  Table  2;  the  values 
of  fractal  dimension  are  derived  from  scanning  electron  micro¬ 
scopic  micrographs  of  the  two  samples. 

The  fractal  parallel  model  (Eq.  (20))  can  be  used  to  pre¬ 
dict  the  upper  bound  of  fceff  of  sample  ‘a’;  the  result  obtained, 
keff  =  1.771  Wm-1  K-1,  is  in  good  agreement  with  the  reference 
value  supplied  by  Toray  Inc.  (1.7  Wm-1  K-1).  With  this  reference 
value,  §  can  be  calculated  as  0.0008  using  Eq.  (23),  which  indicates 
that  the  number  of  perpendicular  channels  in  the  GDL  is  very  small. 
If  the  effect  of  the  perpendicular  channels  is  ignored,  ke  ff  can  be  pre¬ 
dicted  by  the  fractal  parallel  model,  and  an  error  of  4.2%  is  found 
between  the  two  models  for  sample  ‘a’. 


Table  2 

Microstructure  and  material  parameters  of  the  samples  ‘a’  and  ‘b’ 


Sample  ‘a’ 


Sample  ‘b’  Description 


8  x  10-5  m 
3.079  x  10-8  m 
0.78 

1.9  x  10-4  m 
8Wm-'K-' 

0.02624  Wm-1  K_1 
0.25 Wm-1  K  ' 

1.9669 

1.1447 


7  x 10-5  m 
1.487  x  10“8  m 
0.55 

1.9  x  10-4  m 
8Wm->  K-1 

0.02624  Wm-'K'1 
0.25  Wm-1  K_1 

1.9276 

1.1447 


Maximum  pore 
diameter 
Minimum  pore 
diameter 
Porosity 
Thickness  of  the 
gas  diffusion  layer 
Thermal 
conductivity  of 
carbon  fiber 
Thermal 

conductivity  of  gas 

Thermal 

conductivity  of 

PTFE 

Pore  area 

dimension 

Tortuous 

dimension 


fractal  models  and  the  series-parallel  layer  models. 

keff  and  its  upper  bound  of  sample  ‘b’  are  1.7694 Wm-1 1<-1 
and  1.8251  Wm-1 1<-1  calculated  using  Eq.  (29)  with  the  value 
of  §  substituted,  and  Eq.  (27),  respectively.  If  the  effect  of  the 
perpendicular  channels  is  ignored,  kef[  can  be  predicted  by  the 
fractal  parallel  model  (Eq.  (27)),  and  there  is  an  error  of  3.2% 
between  the  two  models. 

The  comparison  between  the  results  of  the  two  samples  shows 
that  the  fceff  value  of  sample  ‘b’  is  greater  than  that  of  sample  ‘a’. 
The  increase  in  keff  is  attributed  to  the  filling  of  PTFE,  which  has  a 
high  thermal  conductivity,  in  the  pore  space  of  carbon  paper. 

Fig.  4  shows  the  prediction  results  of  the  fractal  models  (Eqs. 
(24)  and  (26))  and  the  series-parallel  layer  models  (Eqs.  (5)  and  (6)), 
with  §  =  0.05.  It  can  be  observed  that  the  fractal  model  prediction 
results  are  distributed  between  upper  and  lower  bounds,  indicating 
the  rationality  of  these  models. 

4.2.  Microstructure  parametric  effects  on  effective  thermal 
conductivity 

The  effect  of  the  microstructural  parameters,  the  pore-area 
dimension  Dp,  the  tortuous  dimension  Dt,  and  porosity  e  on  k^ff 
is  studied  using  Eq.  (26)  with  the  data  of  sample  ‘a’. 


Fig.  5.  Effect  of  Dp  on  the  effective  thermal  conductivity  of  carbon  paper. 
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Fig.  6.  Effect  of  Dt  on  the  effective  thermal  conductivity  of  carbon  paper. 


Fig.  5  shows  the  variation  of  k'eff  with  Dp  at  a  fixed  value  of  Dt  =  1.1. 
There  are  two  trends  in  the  variations:  k'eff  decreases  with  increase 
in  Dp  when  ks/kg  >  1,  whereas  it  increases  with  increase  in  Dp  when 
ks/kg  <  1.  It  also  can  be  observed  that  the  effect  of  Dp  on  Id  ff  becomes 
more  significant  when  ks/kg  <  1  than  when  ks/kg  >  1.  On  one  hand, 
the  increase  in  Dp  corresponds  to  the  increase  of  pore  area  of  a 
representative  carbon-paper  unit,  which  leads  to  an  increase  in  the 
proportion  of  the  gas  phase.  On  the  other  hand,  the  relationship  of 
ks/kg  >  1  indicates  that  the  solid  phase  has  a  higher  capacity  of  heat 
transmission  than  the  gas  phase.  With  the  synthesis  of  the  effects  of 
both  Dp  and  ks/kg,  an  important  role  is  detected  for  the  solid  phase 
in  heat  transfer  when  ks/kg  >  1,  and  the  decrease  in  the  proportion 
of  the  solid  phase  with  the  increase  in  Dp  results  in  the  decrease  in 
k'ef{ ;  whereas  the  gas  phase  plays  an  important  role  in  heat  transfer 
when  ks/kg  <  1  and  the  increase  in  proportion  of  the  gas  phase  with 
the  increase  in  Dp  results  in  the  increase  in  k'e{{. 

The  variation  of  k'eff  with  Dt  is  plotted  in  Fig.  6  at  a  fixed  pore-area 
fractal  dimension  of  Dp  =  1.98.  It  can  be  observed  that  fc^ff  decreases 
as  Dt  increases.  The  decrease  in  k'e[[  with  increase  in  Dt  is  due  to 
the  increase  in  thermal  resistance  caused  by  the  longer  distance  of 
heat  transport  through  more  tortuous  channels.  Furthermore,  the 


Fig.  7.  Effect  of  porosity  on  the  effective  thermal  conductivity  of  carbon  paper. 


effect  of  Dt  on  k^ff  becomes  more  significant  when  ks/kg  <  1  than 
when  ks/kg  >  1.  The  reason  is  that  Dt  has  a  significant  effect  on  the 
heat  transfer  of  the  gas  phase,  which  plays  the  leading  role  when 
ks/kg  <  1,  whereas  it  becomes  insignificant  when  ks/kg  >  1. 

Fig.  7  shows  the  effect  of  the  porosity  e  on  k'ef{,  which  is  similar  to 
the  effect  of  Dp.  There  are  two  trends  of  variations  in  this  case  also: 
k^ff  decreases  with  increase  in  e  when  ks/kg  >  1,  whereas  it  increases 
with  increase  in  s  when  ks/kg  <  1. 

5.  Conclusions 

A  novel  fractal  prediction  model  for  the  effective  thermal  con¬ 
ductivity  of  the  GDL  in  a  PEMFC  was  presented  by  considering  the 
real  microstructures  with  the  fractal  theory.  The  prediction  of  the 
fractal  parallel  model  for  carbon  paper,  a  basal  material  of  GDL, 
is  in  good  agreement  with  the  reference  value  supplied  by  Toray 
Inc.  The  prediction  results  from  the  proposed  models  are  also  rea¬ 
sonable  because  they  are  distributed  between  the  upper  and  lower 
bounds.  By  using  the  fractal  series-parallel  model,  the  ratio  of  the 
number  of  perpendicular  channels  to  the  total  number  of  channels 
was  obtained  as  0.0008,  which  indicates  that  few  perpendicular 
channels  exist  in  the  GDL;  therefore,  the  fractal  parallel  model  can 
be  used  for  the  prediction  of  keff  with  a  low  level  of  error.  Moreover, 
the  comparison  between  the  keff  values  of  two  samples  showed  that 
there  isanincreasein  the  fceff  of  GDL  due  to  the  filling  of  PTFE,  which 
has  a  high  thermal  conductivity. 

Finally,  the  microstructural  parametric  effect  on  k^ff  was  stud¬ 
ied  using  the  herein-presented  model  in  dimensionless  formalism. 
It  could  be  concluded  that  k'eff  has  a  positive  correlation  with  the 
effective  porosity  (e)  or  pore-area  fractal  dimension  (Dp)  when 
ks/kg  <  1,  whereas  it  has  a  negative  correlation  with  e  or  Dp  when 
ks/kg  >  1  and  with  tortuous  fractal  dimension  (Dt)  whether  ks/kg  <  1 
or  not. 

The  above  fractal  models  have  the  following  merits:  (1)  The 
equations  deduced  above,  except  Eqs.  (27)-(29),  are  generally 
applicable  for  the  prediction  of  the  effective  thermal  conductivity 
of  two-phase  systems.  (2)  Every  parameter  in  the  model  has  a  clear 
physical  meaning.  (3)  Due  to  the  consideration  of  the  effect  of  pore- 
size  distribution,  this  fractal  model  is  more  realistic  in  comparison 
with  empirical  formulas  and  other  fractal  models.  (4)  By  using 
this  model,  it  is  possible  to  quantitatively  analyze  the  relation¬ 
ship  between  the  microstructure  and  the  macroscopical  material 
parameters  of  the  GDL. 
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